Generalized acceleration theorem for spatiotemporal Bloch waves 
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A representation is put forward for wave functions of quantum particles in periodic lattice po- 
tentials subjected to homogeneous time-periodic forcing, based on an expansion with respect to 
Bloch-like states which embody both the spatial and the temporal periodicity. It is shown that 
there exists a generalization of Bloch's famous acceleration theorem which grows out of this repre- 
sentation, and captures the effect of a weak probe force applied in addition to a strong dressing force. 
Taken together, these elements point at a "dressing and probing" strategy for coherent wave-packet 
manipulation, which could be implemented in present experiments with optical lattices. 

PACS numbers: 72.10.Bg, 67.85.Hj, 42.50.Hz, 03.75.Lm 



I. INTRODUCTION 

The so-called acceleration theorem for wave-packet 
motion in periodic potentials, formulated already in 1928 
by Bloch, ^ has proven to be of outstanding value to solid- 
state physics for understanding the dynamics of Bloch 
electrons within a semiclassical picture.— 12, In its most- 
often used variant, this theorem states that if we con- 
sider an electronic wave packet in a spatially periodic 
lattice, which is centered in k space around some wave 
vector kc-i and if an external electric field E{t) is applied 
under single-band conditions, then this center wave vec- 
tor evolves in time according to hkc(t) — —eE(t), with 
— e being the electronic charge. Perhaps its best-known 
application is the explanation of Bloch oscillations of 
particles exposed to a homogeneous, constant forced"— 
which we recapitulate here in the simplest guise: Take a 
particle in a one-dimensional tight-binding energy band 
E{k) = —{W/2) cos(fca), where W is the band width and 
a denotes the lattice period. Assume that the particle's 
wave packet is centered around fcc(O) initially and sub- 
jected to a homogeneous force of strength F. Then the 
acceleration theorem, now taking the form 



hkcit) = F 



(1) 



tells us kc(t) = kc{0) + Ft/h, so that the packet moves 
through k space at a constant rateii According to another 
classic work by Jones and Zener^ii the particle's group 
velocity Vg{t) in real space is determined, quite gener- 
ally, by the derivative of E(k) with respect to k when 
evaluated at the moving center fcc(t), 



1 dE 

h dfc 



Kit) 



In our case, this relation immediately gives 
Vg{t) = ■^sin(/cc(0)a + WBi) , 



(2) 



(3) 



implying that the particle's response to the constant 
force is an oscillating motion with the Bloch frequency^^ 
lub — Fa/h. This elementary example, to which we 



will come back later in Sec. IIVI strikingly illustrates the 
power of this type of approach. But an obvious restric- 
tion stems from the necessity to remain within the scope 
of the single-band approximation; the above acceleration 
theorem ([Ij is put out of action when several Bloch bands 
are substantially coupled by the external force. Nonethe- 
less, in the present work we demonstrate that there ex- 
ists a generalization of the acceleration theorem which 
can be applied even under conditions of strong inter- 
band transitions. Specifically, we consider situations in 
which a Bloch particle is subjected to a strong oscillat- 
ing force which possibly induces pronounced transitions 
between the unperturbed energy bands. By abandon- 
ing the customary crystal-momentum representation— 
and introducing an alternative Floquet representation in- 
stead, we show that the effect of an additional force then 
is well captured by another acceleration theorem which 
closely mimics the spirit of the original. We obtain two 
major results: The Floquet analog (15^ of Bloch's accel- 
eration theorem ([iJ, and the Floquet analog p2|) of the 
Jones-Zener expression ([2]) for the group velocity. These 
findings are particularly useful for control applications, 
when a strong oscillating field "dresses" the lattice and 
thus significantly alters its band structure, while a sec- 
ond, comparatively weak homogeneous force is employed 
to effectuate controlled population transfer. We first 
outline the formal mathematical arguments in Sees. |TT] 
and mil and then we give two applications of topical in- 
terest, discussing "super" Bloch oscillations in Sec. IIVI 
and coherently controlled interband population transfer 
in Sec. |Vl Although we restrict ourselves here for nota- 
tional simplicity to one-dimensional lattices, our results 
can be carried over to general, higher-dimensional set- 
tings. 



II. THE FLOQUET REPRESENTATION 

We consider a particle of mass m moving in a one- 
dimensional lattice potential V{x) = V{x + a) with spa- 
tial period a under the influence of a homogeneous, time- 
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dependent force F{t), as described by the Hamiltonian 



Zm 



(4) 



Subjecting the particle's wave function ■ip{x, t) to the uni- 
tary transformation 



ip{x,t) =exp( -X / dTF(T) 1 tp{x,t) 



(5) 



the new function ip(x,t) obeys the Schrodinger equation 



d 

ih—tpix^t) = Ho{x,t)^p{x,t) 
with the transformed Hamihonian 



(6) 



Ho{x,t)^^(^p + J^dTFir)^ +V(x). (7) 

Now let us further assume that the force F{t) is periodic 
in time with period T, such that its one-cycle integral 
either vanishes or equals an integer multiple of h times 
the reciprocal lattice wave number 27r/a: 



/ dtF{t) = rxh 
Jq 



27r 
a 



r = 0,±1,±2,... 



(8) 



For example, this is accomplished by a monochromatic 
oscillating force with an additional static bias, 



F{t) = Fr+ Fac COs(cjt) 



(9) 



provided the latter satisfies the condition F^a = rhaj. 
Then the Floquet theorem guarantees that the time- 
dependent Schrodinger equation © admits a complete 
set of spatiotemporal Bloch waves^^"— that is, of solu- 
tions of the form 

ipn,k{x,t) ^ eyip{ikx - ien{k)t/h)un,k{x,t) , (10) 

with spatially and temporally periodic functions 

■Un,k{^^ = "n,fe(2; + a, t) = Un^k{x, t + T) . (11) 

As usual, n is the band index and k a wave number; 
£n{k) thus is the quasienergy dispersion relation for the 
nth band. If r = in Eq. (|8]), the existence of these solu- 
tions is obvious, because then i?o(x,t) = Hq{x + a,t) — 
HQ{x,t + T), so that the wave functions ([ro]) generalize 
the customary Bloch waves ^ for particles in spatially pe- 
riodic lattice potentials by also accounting for the tem- 
poral periodicity of the driving force. When r 7^ 0, so 
that H()[x^t) itself is not periodic in time, spatiotempo- 
ral Bloch waves ([TU)) emerge nonetheless because k is pro- 
jected to the first quasimomentum Brillouin zone, as first 
discussed by Zak.^'^ In any case, the quasienergies e„(fc) 
may depend in a complicated manner on the parameters 
of the driving force, and the wave functions ^n,k{x,t) 



pertaining to a single quasienergy band may be nontriv- 
ial mixtures of several unperturbed energy bands. For 
later use, we observe that their spatial parts 

ifin.kix, t) = exp(ifca;)u„,fc(a;, t) (12) 

obey the quasienergy eigenvalue equation 

d \ 

Ho{x, t) - ih— 1 ipn.k{.X, t) = en{k)ipn.k{x, t) , (13) 

as follows immediately when plugging the solutions 
into the Schrodinger equation ([6]). Throughout, we adopt 
the standard normalization 

27r 

AxLp*^, ^,{x,t)(pn.k{x,t) = Sn,n'6{k - k') . (14) 

An arbitrary wave packet iplxjt) may now be expanded 
with respect to these spatiotemporal Bloch waves, and 
written in the form 



V'(a;,i) = \[^ Jjg'^^ 



t)ipn,k{x,t) , (15) 



with B = [— Tr/a, 7r/a[ denoting the fundamental Bril- 
louin zone. The expansion coefficients gn{k, t) depend on 
the way the system has been prepared and on the way 
the driving force has been turned on, whereas the ba- 
sis functions ipn,k{x,t) and their quasienergies £n{k) are 
given by the eigenvalue equation (jl3p and obviously are 
independent of such details. Clearly, one has 



gn{k,t) = g„(fc, 0)exp( - i£„(fc)t/;i) , 



(16) 



so that the populations \gn{k,t)\'^ remain constant in 
time. This expansion ([T5|). referred to as the Floquet rep- 
resentation of the wave packet, is formally reminiscent of 
its customary crystal- momentum representation, that is, 
of an expansion with respect to the Bloch states of the 
unperturbed potential V{x) which underlies the standard 
acceleration theoremJ^i^i There are, however, substan- 
tial differences which become most clear when consider- 
ing a wave packet occupying a single quasienergy band. 



il^ix,t) 



dkg{k,t)(pk{x,t) ; 



(17) 



here and in the following, we omit the band index n for 
ease of notation. Now this wave packet (fTT]) may de- 
scribe, for instance, the dynamics in a situation where 
two unperturbed energy bands are resonantly coupled 
by the driving force F{t); consequently, in a crystal- 
momentum representation one would have to account 
for Rabi-type oscillations between these two bands by 
coefficients which quantify the oscillating band popula- 
tions. In the Floquet respresentation, on the other hand, 
the Rabi oscillations are already incorporated into the 
basis states (ITUl) . so that one merely encounters single 
quasienergy band dynamics, with the remaining time 
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evolution of g{k,t) simply given by Eq. (IT6l) . Thus, 
although the external force effectuates transitions be- 
tween the unperturbed Bloch bands, there are no inter- 
quasienergy band transitions; \g{k,t)\'^ remains constant 
in time. Second, even in a situation where F{t) does not 
couple different energy bands, the wave packet's center 
kc {t) evolves according to the standard acceleration the- 
orem hkc = F in the crystal-momentum representation, 
whereas in the Floquet representation the moment 



(fc> 



dkk\g{k,t)\' 



(18) 



obviously stays constant in time. In short, an expansion 
of the wave packet with respect to the spatiotemporal 
Bloch waves ([T0| implies constant coefficients, and hence 
constant occupation probabilities, if the external force 
F{t) adheres to the specification ([8|). This formal shift of 
the dynamics from the occupation numbers to the basis 
states which is implied by the Floquet representation now 
allows for a clear and physically transparent description 
of the additional effects which emerge when the external 
force does not obey Eq. ([S]); these effects are captured 
by the generalized acceleration theorem exposed in the 
following. 



III. THE FLOQUET ACCELERATION 
THEOREM 



We take a wave packet occupying a single quasienergy 
band and stipulate that in addition to the possibly strong 
driving force F(t) there is a second homogeneous force 
Fp(t) which we denote as the probe force; this is assumed 
to be sufficiently weak so that it does not introduce tran- 
sitions among different quasienergy bands. To be precise, 
the total Hamiltonian now reads 



H{x,t) 



2m 



+ V{x)- F{t)x- Fp{t)x , 



(19) 



where the time-periodic force F{t) is resonant in the sense 
of Eq. ([U and thus creates a basis of spatiotemporal 
Bloch waves (|10|) . whereas the probe force Fp{t) also is 
spatially homogeneous, but not necessarily periodic in 
time. After performing the unitary transformation ([5]), 
we obtain the Hamiltonian in the form 



H{x,t)=Ho{x,t)-Fp{t)x , 



(20) 



with Ho(x^t) given by Eq. ([7]). Moreover, we start from 
an initial wave packet of the form ([T7|) . Because of the 
additional probe force Fp{t), the time evolution of g{k, t) 
is no longer given by Eq. ()16|) : the aim now is to find an 
effective Hamiltonian H which governs the resulting dy- 
namics oi g{k,t), under the proposition that this remains 
restricted to the single, initially occupied quasienergy 
band. 

Exploiting the normalization (fT4|) . we have 



dxipl{x,t)ip{x,t) 



(21) 



This gives 



a 







Ho - ift| 





—Fp I dxLplx^p 



(22) 



having suppressed the arguments x and t for better legi- 
bility; all integrals here are taken over the entire lattice. 
In the first term on the right-hand side of this equation 
we exploit the quasienergy eigenvalue equation, Eq. (|13l) , 
yielding e{k)g{k, t). For rewriting the second term we use 



(23) 



which is obtained by taking the derivative of the complex 
conjugate to Eq. ([T2|) with respect to k, and leads to 



^ j dxifilx^p = idkg " iy'^ J dxe~ 
idkg - i{dkUk\uk)g 



ikx 



(24) 



For making the final step, we have resubstituted the ex- 
pression ()17|) for tJj and have made use of the identity 



' dxe'^'''-''^^uk'dkul = —Sik - k'){dkUk\uk) , (25) 
a 

with the scalar product 



{dkUk\uk) ^ / dxuk{x,t)dkul{x,t) (26) 
Jo 

being given by an integral over a single lattice period. 
Note that 

{uk\uk) = 1 , (27) 
as an immediate consequence of Eq. (jl4p , which implies 

{dkUk\uk) + {dkUk\uk)* ^ Q , (28) 

so that {dkUk\uk) is purely imaginary. Collecting all the 
pieces, we obtain the desired evolution equation 



d 

^fi-g^g{k,t) = ng{k,t) , 



(29) 



with the effective Hamiltonian for the Floquet represen- 
tation, 



v. = e{k) - iFpdk - Fp lui{dkUk\uk) 



(30) 



From this expression we deduce the generalized acceler- 
ation theorem, that is, the acceleration theorem for the 
Floquet representation: Since the moment (|18p obeys the 
the equation 



dt 



(fc) 



1 



[H,fc]) 



(31) 
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and the commutator appearing here on the right-hand 
side is easily evaluated, i[H,k] = Fp, we are directly led 
to 



On the other hand, exploiting the operator identity 



h-{k)it) ^ Fpit) . 



(32) 



This is the central result of the present work; its analogy 
to the standard acceleration theorem ([T]) for the crystal- 
momentum representation is evident. Observe that there 
is an intuitively clear reason for the appearance of the 
term proportional to (9fcWfc|ufc) in the effective Hamil- 
tonian ([30|): The twofold periodic parts Uk{x,t) of the 
spatiotemporal Bloch waves are obtained by solving the 
eigenvalue equation (jl3p . This is done for each wave 
number k separately, so that one is free to bestow upon 
each eigensolution an arbitrary phase factor exp(i0(fc)). 
On the other hand, the evolution equation ([29]) for the 
wave function g{k,t) in the Floquet representation nat- 
urally establishes a "connection" between those differ- 
ent eigensolutions^i^ and therefore requires information 
about the gauge function 0{k); this is provided by the ex- 
pression {dkUk\uk)- Note further that when multiplying 
Eq. (1^ by g*{k,t) and subtracting the complex conju- 
gate of the resulting equation, this piece drops out, and 
one is left with 



d Fpjt) d 
dt h dk 



(33) 



Thus, \g{k,t)\'^ does not depend on k and t separately, 
but rather on the combination k — J^dr i^p(r)//i, so that 
the distribution g(k,t) moves through the Floquet k 
space without change of shape, again in precise anal- 
ogy to the classic behavior.^ But we reemphasize that 
this seemingly simple dynamics might be unrecognizable 
in the usual crystal-momentum representation, because 
the system might undergo violent transitions between dif- 
ferent energy bands when monitored in a basis of time- 
independent Bloch waves. 

As the introductory example has shown, the stan- 
dard (crystal-momentum) acceleration theorem develops 
its main power in combination with the Jones-Zener ex- 
pression (21) for the wave packet's group velocity in real 
space, and the question naturally arises whether there 
exists a similar connection in the Floquet representation. 
Obviously, one can establish a relation corresponding to 
Eq. ([2]) by applying a stationary-phase argument to the 
expansion (jlSp . but here we follow an alternative line of 
reasoning which may be found particularly enlightening. 
Considering a well-localized wave packet ip{x,t) in the 
original frame of reference to which the Hamiltonian op- 
erators (IH) and ([TO)) pertain, that packet's group velocity 
is given by 



^sit) = 3-(V'(a;,i)k|V'(a:,t)) 



dt 
1 

m 



(35) 



the eigenvalue equation ()13p transforms into the even 
more basic eigenvalue equation 



d 

Hk{x, t) - ih— ) Un,k{x, t) = en{k)un,k{x, t) (36) 



for the periodic core pieces Un,k(x,t) of the spatiotem- 
poral Bloch waves (fTO)) . invoking the parametrically k- 
dependent operator 

Hk{x,t)^^(^p + nk + J^dTF{T)^ +Vix). (37) 

This eigenvalue problem can efhciently be implemented 
for numerical calculations)^ It also manifestly contains 
the origin of the condition ([SJ imposed on the oscillating 
force F{t), since k is reduced to fall within B. Most im- 
portantly, this eigenvalue problem ([55)) poses itself in an 
extended Hilbert space made up of functions Un^kix,t) 
which are periodic in both space and time, in accordance 
with Eq. ([TT)) . Consequently, "time" has to be regarded 
as a coordinate in this extended Hilbert space and there- 
fore needs to be integrated over when forming a scalar 
product, just like any spatial coordinate. Thus, the nat- 
ural scalar product in this extended Hilbert space is given 
by27 



(38) 



with ( • I • ) denoting the standard scalar product in the 
original, physical Hilbert space, as already employed in 
Eqs. ([26| and ([27| . It follows that the quasienergies £„(fc) 
can be written as diagonal elements of the matrix of the 
quasienergy operator. 



£„(fc) = {{UnM\Hk - ifldt\Un,k)) 



(39) 



{^{x,t)\p\^Pix,t)) . 



(34) 



inviting us to make use of an analog of the Hellmann- 
Feynman theoremj^l 

-|^£„(fc) = {{Un,k\'^^\Un,k)) 

dk dk 

h /■* 

= —{{un.k\p + fi-k+ / dTF(r)|u„,fc)) 
m Jo 

- -((^„,fcb|V^n,fc)) . (40) 

m 

In the final step we have undone the shift (1551) : the wave 
functions ipn,kix,t) then denote the functions which are 
obtained from the spatiotemporal Bloch waves ([TO)) by in- 
verting the transformation ([5]). Comparison of Eqs. ([M)) 
and ([40)) . keeping in mind the definition ([38)) . now yields 

the desired relation: Supposing that ip{x,t) = ipn,koi^'^) 
were made up from a single spatiotemporal Bloch wave 
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labeled by n and ko, say, one would obtain the formal 
identity 



1 

- J dtv^it) = 



1 de^ 
h dk 



(41) 



But this is not what we want, because an individual 
spatiotemporal Bloch wave is uniformly extended over 
the lattice and thus does not correspond to a "group" 
which propagates in space. Rather, we require a wave 
packet (fT7|) which is reasonably well centered in the Flo- 
quet k space, with a center (fc) given by Eq. (fT8|). Then 
we have 



1 de^ 
h dk 



(42) 



(fe) 



to good accuracy, so that the cycle-averaged group veloc- 
ity of the Floquet wave packet is given by the derivative of 
its quasienergy dispersion relation, evaluated at its center 
(fc). Again, this Floquet relation (|42|) closely mimics its 
historic crystal-momentum antecessor, given by Eq. 
In contrast to the equation of motion ((32)) for (fc) itself, 
which holds exactly within a single quasienergy band set- 
ting, this relation is an approximation which holds 
the better, the narrower the packet's Floquet k space 
distribution. Although it seems self-evident, it might 
be worthwhile to stress that the argument required to 
evaluate the derivative (H^ is Floquet (fc), not crystal 
momentum kc- 



IV. SUPER BLOCH OSCILLATIONS 

The phenomenon termed "super" Bloch oscillations^^ 
arises when a Bloch particle is subjected to both a static 
(dc) and an oscillating (ac) force, such that an integer 
multiple of the ac frequency is only slightly detuned from 
the Bloch frequency associated with the dc component of 
the force 1^"— Although the effect itself appears almost 
trivial from the mathematical point of view, we nonethe- 
less dwell on this at some length, because it provides a 
particularly instructive example for juxtaposing the fa- 
miliar crystal-momentum representation to the Floquet 
representation introduced in Sec. Ulland for demonstrat- 
ing in detail how they match. To be definite, we consider 
the total force to be of the form 



F{t) = e(t - to) [Fdc + Fac cos{ujt)] 



(43) 



where Q{t) denotes the Heaviside function, so that both 
the dc and the ac component of the force are turned on 
instantaneously and simultaneously at to; that moment 
to thus determines the relative phase between the Bloch 
oscillations caused by the dc component and the driving 
oscillations of the ac component. 

The basic assumptions now are that (i) we are given an 
initial wave packet which occupies a single energy band, 
being centered around fcc(io) at the moment t = tg in 



the crystal-momentum representation, and that (ii) in- 
ter band transitions remain negligible for t > to, despite 
the action of the force F{t). We then encounter single- 
band dynamics which are fully covered by the "old" ac- 
celeration theorem hkc{t) ~ F{t), giving 



kc{t) 



kc{to 
1 



(44) 



Fdc{t — to) H — sin(a;t) — sai{Loto) 



for t > to- As an archetypal example we now take a tight- 
binding cosine energy dispersion relation for the band 
considered. 



w 



cos(fca) , 



(45) 



parametrized as in the Introduction. Utilizing Eq. ([2l). 
one then finds the packet's group velocity: 



I dE 

- n dk 



feo(t) 



— sm(fcc(t)a) 



(46) 



This expression describes super Bloch oscillations if we 
assume further that the dc component of the force is 
almost resonant in the sense of Eq. We therefore 

decompose this component according to 



i^dc =Fr + SF 



(47) 



where FrU = rhu with some nonzero integer r as previ- 
ously in Eq. ([9]), so that rw equals the Bloch frequency 
wb = Fr-a/h, while 5F is quite small compared to F^. We 
then have 



(48) 



with frequency detuning 5uj = SFa/h, so that the group 
velocity ([46]) takes the form 

Vg{t) = sm(rujt + Sujt + K sm{ujt) + $^ , (49) 
having introduced the scaled driving amplitude 



and a constant phase 

•I* = kc(to)a — (j'o; + Suj)to — K sm{ujto) 



(51) 



which accounts for the initial conditions. Because 6uj <C 
cj according to our specifications, the contribution Scot 
to the argument of Vg does not vary appreciably during 
one single cycle T = 2tt/uj of the ac component. Thus, 
when averaging the instantaneous group velocity over one 
such cycle, this "slow" time dependence may be ignored, 
meaning that 6u! t may be considered as constant when 
taking the averagei^ Invoking the Jacobi- Anger inden- 
tity in the guise 



(52) 



6 



where Ji{K) denote the Bessel functions of the first kind, 
one immediately obtains 



1 

7f dtv^it) 



(53) 



According to the above reasoning, here the "fast" time 
dependence is integrated out, but the slow dependence on 
6uj t remains.— Integrating, this yields the cycle-averaged 
drift motion of the packet, that is, its position Xg{t) with- 
out the fast ac quiver, 

W 

x,{t) = -{-iyj,{K)—cosi6L. t + $) , (54) 

with a suitably chosen origin of the x axis. This result 
finally clarifies what is "super" with these dynamics: Be- 
cause the residual force SF is quite small, the amplitude 
of this oscillation (|54|) can be fairly large; indeed, in a 
corresponding experiment with weakly interacting Bose- 
Einstein condensates in driven optical lattices Haller et 
al. have observed giant center-of-mass oscillations with 
displacements across hundreds of lattice sites»^ 

As far as the phenomenon itself is concerned there is 
nothing more to add; because one requires single Bloch- 
band dynamics right from the outset, a Floquet treat- 
ment is not necessary. Nevertheless the Floquet approach 
is of its own intrinsic value even here, since it provides a 
somewhat different view which, in contrast to the above 
crystal-momentum calculation, is capable of some gener- 
alization. 

The Floquet analysis starts from the spatiotemporal 
Bloch waves and their quasienergies. In a single-band 
setting with an external homogeneous force, these are 
exceptionally easy to obtain: Writing the Bloch waves of 
the undriven lattice in the form 



^k{x) 



' wi{x)e' 



Ma 



(55) 



where wi{x) denotes a Wannier function localized around 
the ith lattice site^ the so-called Houston functions^^ 

V^fe(x,i) = ^w;,(a;)e'«'=(*)^''exp('-i / drS(gfc(r)) 

(56) 

are solutions to the time-dependent Schrodinger equation 
in the original frame, for arbitrary F{t), provided the 
"moving wave numbers" qk{t) are given by 



qk{t)^k+^j^ATF{T) 



(57) 



always assuming the viability of the single-band approx- 
imation.'^* Taking a force of the particular form ^ with 
exactly resonant F^ obeying Fj-a = rhhj, we have 



IkV-) — K -\ h -r — sm(a;n 

a nuj 



(58) 



This implies that both the exponentials exp(iqk{t)£aj 
and E{qk{t)) are T periodic in time, with T = 2tt/uj, 
whereas the integral over E[qk{t)) is not, because the 
Fourier expansion of E(^qk(t)^ contains a zero mode, so 
that its integral contains a linearly growing contribution. 
But this observation reveals that the "accelerated Bloch 
waves" (j56p with resonant time-periodic forcing Q are 
precisely the required spatiotemporal Bloch waves in the 
original frame, with their quasienergies being determined 
by the zero mode: 



e{k) 



T 



dtE{qk{t)) 



W 



= -(-l)'^J,.(i^)y cos(fca) 



(59) 



The remarkable fact that the quasienergy bands collapse, 
i.e., become flat when K is such that Jr{K) = 0, indi- 
cates that an oscillating force can effectively shut down 
the tunneling contact between neighboring wells; this 
"coherent destruction of tunneling" is a generic feature of 
driven single-band systemsi^"— A bit of reflection then 
shows that the core pieces Uk{x,t) of the spatiotempo- 
ral Bloch waves, that is, the solutions to the eigenvalue 
equation (j36p . are given by 



Uk 



(x) = ^^«,(x)e'"" 



X expl dr [E{qk{T)) ~ e{k)yj . (60) 



Although this has not been particularly emphasized, the 
above construction makes sure that any spatiotemporal 
Bloch wave (j56l) is labeled by the same wave number 
k as the ordinary Bloch wave to which it reduces when 
the external force vanishesi^ Otherwise, there is nothing 
particular about the choice t = for the lower bound of 
integration in Eq. ([57|) for q^ {t) : In contrast to Eq. (|43|) , 
where t = to has been singled out as the moment when 
the force is turned on, and which thus designates an 
initial- value problem for a particular wave packet, the 
solution of the eigenvalue problem ([55)1 for the entire 
spatiotemporal Bloch basis requires a force F{t) which 
is perfectly periodic in time; the resulting expression for 
qk{t) thus holds for both t > and t < 0. Also note that 
it would be meaningless to include some additional con- 
stant phase into the argument of the ac component of the 
force Because this expression holds for all times t, 
such a phase would merely amount to a shift of the origin 
of the time coordinate and thus is as irrelevant for the cal- 
culation of the quasienergy dispersion relation as would 
be a shift of the origin of the spatial coordinate system 
for the calculation of a crystal's energy band structure. 

Knowing the quasienergy dispersion relation (|59p . the 
machinery established in Sec. Illll can be set to work: Ac- 
cording to Eq. (jHI, the cycle-averaged group velocity of 
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a Floquet wave packet p7|) is given by 



1 de 
h dk 



(k) 



(61) 



If we now turn back to the specific forcing (|43l) , and thus 
consider exactly the same initial- value problem as in the 
previous crystal-momentum exercise, we can make op- 
erational use of the decomposition (|T7)) of the dc force: 
Its resonant part has already been incorporated into 
the spatiotemporal Bloch waves ([56)) . which means that 
it has already been accounted for in "dressing" the lat- 
tice and changing its original energy dispersion E{k) to 
the quasienergy dispersion e(fc). Therefore, it is only the 
small residual part SF which enters into the equation of 
motion for (fc), that is, into the generalized acceleration 
theorem p2l) : this part SF thus constitutes a particu- 
lar, time- independent example of a probe force Fp{t) as 
considered in Sec. IIIII We now have 



givmg 



^ d , , , , , hSuj 
dr a 



{k)it) = {k)it„) + —it-h). 



(62) 



(63) 



All that remains to be done now is to express the ini- 
tial Floquet center {k){to) in terms of the initial wave 
packet's center kc(to), which had been specified in the 
crystal-momentum representation. But this is an easy 
task, comparing the original Bloch waves ([55]) to their 
spatiotemporal descendents (|56|) : At the moment Iq when 
the force (|^ is turned on, fcc(io) coincides with qk{to) 
for one particular k] this evidently is the desired (fc)(to)- 
The equality identifying {k){to) thus is 



kc{to) = 9(fc)(to)(^o) , 
which, written out in full detail, reads 



(64) 



fcc(io) = (fc)(io) + + sin(c.io) • (65) 
a nuj 

Using this to eliminate {k){to) from Eq. (|63| . we arrive 
at 



bit mind-boggling; if this is not done with sufficient care, 
one might overlook a contribution to ^."^^ But if respected 
properly, the mathematical structure of the Floquet pic- 
ture unerringly leads to the correct answer. 

If one strips the above reasoning to the bare essentials, 
that is, if one starts from the quasienergy dispersion re- 
lation (15^ . takes its derivative to obtain the formal ex- 
pression (|6T|) for the cycle-averaged group velocity, and 
then inserts the solution to the equation of motion ([5^ 
dictated by the generalized acceleration theorem in or- 
der to compute the group velocity of the wave packet 
actually considered, one sees that this procedure exactly 
parallels the explanation of the usual Bloch oscillations, 
as reviewed in the Introduction. Thus, super Bloch oscil- 
lations may be seen as ordinary Bloch oscillations arising 
in response to a weak probe force 6F, but occurring in a 
spatiotemporal lattice, as created by dressing the original 
lattice through application of the strong force 

One might finally wish to get away from the particular, 
instantaneous onset of the forcing assumed in Eq. (|43l) : 
The dc and the ac component might not be switched on 
simultaneously, or not abruptly, possibly involving two 
different turn-on functions for the two components. In 
any case, at some moment to the final amplitudes will 
have been reached, so that the previous analysis goes 
through unaltered for t > Iq, it one only interprets fcc(io) 
correctly: This would no longer indicate the crystal- 
momentum wave number around which the initial wave 
packet had been prepared, but rather that to which the 
latter had been shifted during the turn-on phase. Ex- 
pressed differently, the phase $ in Eqs. (|53|) and (|54|) de- 
pends significantly on the precise turn-on protocol: Not 
surprisingly, the way the external force has been turned 
on in the past crucially affects the coherent wave-packet 
motion after the turn-on is over. 

Aside from its aesthetic value, the Floquet picture of- 
fers at least one further benefit: Bloch oscillations in 
dressed lattices may also occur under conditions such 
that the quasienergy bands are mixtures of several un- 
perturbed energy bands, disabling a crystal-momentum 
treatment. A Floquet analysis, on the other hand, would 
merely require one to replace the single-band quasiener- 
gies ([5^ by the actual ones and then again invoke the 
generalized acceleration theorem ((32|) . similar to the ex- 
amples worked out in the next section. 



{k){t)a = kc{tQ)a — rujto — K sm{iotQ) + Suj(t — to) 

= 5u}t + <i> (66) 



COHERENT CONTROL OF INTERBAND 
POPULATION TRANSFER 



with precisely the same phase $ as already defined in 
Eq. (fSTj) . Inserting this argument (|66| into the cycle- 
averaged group velocity ()61|) . and comparing with the 
previous expression ([55)1 . one confirms that the result 
of the Floquet analysis fully coincides with that of the 
more customary crystal-momentum calculation. The ne- 
cessity to painstakingly distinguish between crystal mo- 
mentum kc and Floquet (fc) at all stages may appear a 



A field of major current interest in which the Floquet 
picture may find possibly groundbreaking applications 
concerns ultracold atoms, or weakly interacting Bose- 
Einstein condensates, in time-periodically driven opti- 
cal lattices.— i^i^i^i^"— As opposed to ordinary crys- 
talline matter exposed to high-power laser fields, such 
systems offer the advantage that one can apply even non- 
perturbatively strong driving forces without inducing un- 
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wanted inhomogeneities, as caused by polarization effects 
or domain formationJS, The issue at stake here is not 
merely redoing well-known condensed-matter physics in 
another setting, and thus selling old wine in new skins, 
but rather finding genuinely new ways of coherently con- 
trolling mesoscopic matter waves, such that target states 
are created which have not been accessible before, and 
are manipulated according to some prescribed protocol. 
Here we point out that the generalized acceleration the- 
orem ([5^ may be a valuable tool in this quest. 

A standard one-dimensional (ID) optical lattice is de- 
scribed by a cosine potential 



V{x) — — coa{2ki^x) 



(67) 



where kj^ is the wave number of the two counterpropagat- 
ing laser beams generating the lattice.— Its depth Vq is 
measured in multiples of the single-photon recoil energy 



2m 



(68) 



For orientation, if one traps *^Rb atoms in a lattice with 
/cl corresponding to the wavelength A — 842 nm, as in 
a recent experiment by Zenesini et al.^ one finds E,- — 
1.34 X 10~^^ eV; typical optical lattices are a few recoil 
energies deep. 

Figure [T] shows quasienergy spectra for such a ID co- 
sine lattice (1571) with depth Vo/E^ — 5.7 under pure ac 
forcing, that is, for F{t) — F^c cos{ujt) not containing 
a dc component, with driving frequency uj = 3.71 Ey/h. 
Under the laboratory conditions specified above (^^Rb at 
A = 842 nm), this corresponds to u}/{2tt) — 12 kHz. Fig- 
ure [TJa) results when the scaled driving amplitude (|50)) 
is set to zero; this subfigure therefore is obtained by pro- 
jecting the lowest three unperturbed energy bands to 
the fundamental quasienergy Brillouin zone, which ex- 
tends from e = —hbj /2 to e = +fiaj/2 on the ordinate. 
Figure [Ijb) displays the quasienergy band structure for 
the moderate driving strength K — 0.5; here avoided 
crossings show up which generally indicate multiphoton- 
like resonances. Figure[ljc) then reveals pronounced ac 
Stark shifts (that is, shifts of the quasienergies against 
the zone-projected original energies) for K — 3.0, corre- 
sponding to truly strong forcing. 

We now turn from the quasienergy spectrum to an ex- 
emplary initial-value problem: At < = we prepare an 
initial wave packet p7)) in the lowest Bloch band n = 1 
with a Gaussian momentum distribution. 



gi(fc,0) = (y¥A/c)-i/2exp 



2{Akf 



(69) 



centered around fcc(0)/fcL = with width Afc/fcL = 0.1, 
and subject it to a pulse. 




FIG. 1; (a) Quasienergy spectrum of an ac-driven ID opti- 
cal lattice (|67|l with depth Vo/Er — 5.7, scaled driving fre- 
quency fkd/Er — 3.71, and scaled driving amplitude K — 0. 
This figure is obtained by projecting the lowest three energy 
bands of the undriven lattice to the first quasienergy Brillouin 
zone, ranging from e/(Jiw) = —1/2 to e/{ttw) — +1/2. (b) 
Quasienergy band structure for K = 0.5. Here the ac Stark 
shifts still are comparatively weak, but the time-periodic forc- 
ing introduces pronounced avoided crossings among the "low- 
est" three bands, (c) Quasienergy band structure for K — 3.0 
for the "lowest" five bands, revealing substantial ac Stark 
shifts. 



starting at t = and ending at t = Tp, endowed with a 
smooth, squared-sine envelope function: 



sit) 



sm 



^) , 0<t<Tp. (71) 



We again set uj = 3.71 E,-/h, as in Fig.[T] adjust the pulse 
length to 50 cycles; Tp = 50 x 27r/w, and fix the maximum 



driving amplitude -Fmax such that 



^a/{hiu) 



F{t) = F^^^sit) sin(wi) 



(70) 



3.0, corresponding to the conditions reached in Fig.[IJc). 
We then monitor the resulting wave-packet dynamics 
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25 50 

t/T 

FIG. 2: (Color online) Time evolution of an initial wave 
packet (|69p under the action of a driving pulse (|70p with the 
smooth squared-sine envelope \71\ . with maximum scaled am- 
plitude A'niax = 3.0, scaled driving frequency hw/E^ = 3.71, 
and pulse length Tp = SOT, where T — 2n/ijj is the dura- 
tion of a single cycle. The main frame shows the occupation 
probabilities of the original Bloch energy bands during the 
pulse. Apart from the initially populated lowest band n = 1 
(jagged line at the top), both bands n = 2 and n — 3 become 
significantly excited during the pulse (jagged lines at the bot- 
tom), the band n = 3 even to a higher extent than n — 2 
at maximum driving strength. In contrast, when monitoring 
the same dynamics within the bases provided by the instan- 
taneous spatiotemporal Bloch waves, only the corresponding 
Floquet band n = 1 is appreciably occupied, as shown by the 
horizontal line at the top and magnified in the inset. Observe 
the scale of the inset's ordinate! 



both in the basis of the unperturbed energy bands and 
in the bases provided by the instantaneous spatiotempo- 
ral Bloch waves, that is, in the family of Floquet bases 
which are obtained when the driving amplitude is kept 
fixed at any value Fac = i^maxs(io) reached during the 
pulse. Figure [5] displays the results: The jagged lines in 
the main frame show the occupation probabilities of the 
lowest three unperturbed Bloch bands n = 1,2, and 3; 
in the middle of the pulse the band n — 3 contains even 
more population than the band n — 2. On the other 
hand, the horizontal line at the top depicts the occu- 
pation of the instantaneous Floquet band emerging from 
the lowest Bloch band: This Floquet band contains prac- 
tically all the population during the entire pulse, which 
means that the wave function adjusts itself adiabatically 
to the changing morphology of its quasienergy band,'^^ as 
previously sketched in Fig. [TJ when the driving amplitude 
^maxs(i) is first increased and then decreased back to 
zero. To quantify the precise degree of adiabatic follow- 
ing, the inset in Fig. [5] shows the variation of the Floquet 
band population on a much finer scale. Observe that the 
final adiabaticity defect is on the order of merely 0.1%, 
even though the driving amplitude reaches its fairly high 
maximum strength within no more than 25 cycles. 

With respect to the concepts developed in Sec. HIl 
Fig. [2] strikingly demonstrates the advantages of the Flo- 
quet picture over the traditional crystal-momentum rep- 



resentation for the situation considered. If there were 
an additional probe force, its effect would have to be 
tediously disentangled from the fast oscillations of the 
Bloch band populations. When the same dynamics are 
seen from the Floquet viewpoint, essentially "nothing" 
happens, because practically all inter-Bloch-band transi- 
tions are already accounted for by continuously adapt- 
ing the Floquet basis, so that the action of a probe 
force would stand out most clearly. Although, of course, 
the crystal-momentum representation is mathematically 
equivalent to the Floquet picture, there is no question 
which one is preferable here. Note also that Fig. [5] an- 
swers one further pertinent question: How do we prepare 
a wave packet which occupies merely a single quasienergy 
band, although it is undergoing rapid transitions between 
several Bloch bands at the same time? The recipe for 
achieving this is simple: Start with a wave packet occu- 
pying a single Bloch band and switch on the driving force 
smoothly thereby enabling adiabatic following. 

At this point an important issue needs to be stressed: 
The concept of adiabatic following, or parallel trans- 
port in a differential-geometric language, usually is ap- 
plied to energy eigenstates j^^i^^ in the context of opti- 
cal lattices this has been exploited, e.g., by Frataloc- 
chi and Assanto for studying nonlinear adiabatic evolu- 
tion and emission of coherent Bloch waves4^ In contrast, 
here we consider adiabatic following of explicitly time- 
dependent quasienergy eigenstates, that is, of solutions 
to the quasienergy eigenvalue equation (|13l) ; this is what 
allows us to separate the fast, oscillating time depen- 
dence of the driving force from the slow, parametric time 
dependence of its envelope. 

Having learned these lessons, we now set the gen- 
eralized acceleration theorem (1321) to work. Suppose 
that we are prompted to empty the ground-state energy 
band. Starting again from an initial wave packet (j69l) . 
we then may proceed as follows: First we smoothly turn 
on an ac force which dresses the lattice, creating avoided 
quasienergy crossings initially not "seen" by the adia- 
batically following packet. For instance, we may wish 
to utilize the avoided crossings showing up in Fig. (rfb). 
To this end, we again take an ac force with frequency 
uj = 3.71 E^/h and fix its scaled driving amplitude at the 
plateau value K = 0.5. This dressing force is switched 
on during 25 cycles with half a squared-sine envelope, 
maintained at maximum amplitude for 50 further cy- 
cles, and switched off again for another 25 cycles, as 
sketched in Fig. EJa). If this were all we did, the wave 
packet would simply undergo adiabatic evolution and fi- 
nally restore its initial condition, as previously observed 
in Fig. [2] Instead, once the maximum dressing ampli- 
tude has been reached, we now apply an additional weak 
probe force Fp{t) in order to exploit Eq. ([5^ for mov- 
ing the packet away from the Brillouin zone center, driv- 
ing it over the avoided crossings that have opened up in 
Fig. mb). This probe force is implemented in the form 
of two smooth, squared-sine shaped dc pulses, one acting 
during the plateau of the dressing pulse, the other act- 
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FIG. 3: (Color online) (a) Protocol for achieving almost com- 
plete interband population transfer in a dressed optical lattice 
by means of a weak probe force. The dashed line is the scaled 
envelope of the dressing ac force. The solid line is the neg- 
ative, scaled probe force, amplified by a factor of 10. (b) 
Resulting wave-packet dynamics in the Floquet representa- 
tion, shown as a contour plot of Under the first 
action of the probe force the wave packet is shifted in accor- 
dance with the generalized acceleration theorem (|32p , until it 
undergoes Zener-type transitions to other quasienergy bands 
at the avoided crossings visible in Fig. [IJb). After the dress- 
ing force is switched off, the second, reversed action of the 
probe force shifts the remaining part of the packet back to 
the Brillouin zone center, (c) Comparison of the initial wave 
packet (dashed line) with the part of the wave function that 
remains in the lowest energy band at the end of the process 
(solid line). 



ing with reversed sign after the dressing pulse is over, as 
drawn in Fig. Efa) . The maximum strength of the probe 
force here is only 2.5% of that of the dressing force; for 
better visibility, the probe force is magnified in Fig. 
by a factor of 10. 

It is now almost obvious how to describe the response 
of the wave packet within the Floquet picture: The ini- 
tial state (IMj) first is adiabatically shifted into a sin- 
gle quasienergy-band packet during the turn-on of the 
dressing force. In contrast to a crystal-momentum repre- 



sentation, all dressing-induced fast oscillations are taken 
out of the dynamics of gi{k,t) in the Floquet represen- 
tation, as shown in Fig. EJb). When the first probe 
pulse acts at constant dressing amplitude, it forces the 
wave packet over the avoided crossing seen in Fig. [I{b), 
so that the packet undergoes Zener-type transitions to 
"higher" quasienergy bands, '^^'^^ splitting into individ- 
ual subpackets associated with the different quasienergy 
bands involved. When the dressing force is switched 
off, each of these subpackets moves adiabatically on its 
own quasienergy surface, finally reaching the continu- 
ously connected Bloch bands. The second, reversed 
probe pulse, applied after the dressing pulse is over, then 
acts in accordance with Bloch's original acceleration the- 
orem ([T]), shifting the various subpackets back to the Bril- 
louin zone center. In the scenario displayed in Fig. |31 the 
lowest band is almost entirely depopulated by the probe- 
induced Zener transitions, so that only a marginal frac- 
tion of the initial packet returns, as depicted in Fig.[3l^c). 
Thus, the main part of the initial packet has been placed 
in higher Bloch bands, as intended. We have also checked 
by explicit calculation that without the comparatively 
weak probe pulses the returning wave packet would be 
almost identical to the initial one. 

The above example of our "dressing and probing" 
strategy immediately lends itself to a host of further 
modifications and extensions. To give but one further 
instance, if the probe pulse is still weaker, such that 
the wave packet does not pass over the avoided-crossing 
regime, but rather stops there, the Zener transitions are 
incomplete, so that a signifcant part of the initial state 
is recovered when the process is over. This is elaborated 
in Fig. 21 with the same dressing force as above, but now 
the maximum strength of the probe force amounts to only 
1.7% of that of the dressing force. The final subpacket 
still occupying the lowest Bloch band then is no longer 
centered around fc/^L — 0, implying that this subpacket 
will move over the lattice. In a sense, the left wing of the 
initial wave packet has been cut out, so that Fig. |4]may 
be regarded as a particular paradigm of "wave-packet 



surgery. 



"44 



VI. CONCLUSIONS 

Summarizing our line of reasoning, we have introduced 
in Sec. HI] a representation of wave packets of quantum 
particles in spatially periodic lattices subjected to ho- 
mogeneous, time-periodic forcing which is based on an 
expansion with respect to spatiotemporal Bloch waves 
and reduces to the standard crystal-momentum repre- 
sentation when the forcing is turned off. It embodies 
forcing-induced oscillations into the basis, so that only 
the actually relevant dynamics remain to be dealt with. 
Within this Floquet representation one encounters many 
features already familiar from solid-state physics in time- 
independent lattice potentials, but here their scope is dif- 
ferent. As a prominent example, the generalized acceler- 
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FIG. 4: (Color online) (a) Protocol for achieving partial inter- 
band population transfer in a dressed optical lattice by means 
of a weak probe force. The dashed line is the scaled envelope 
of the dressing ac force, which is the same as in Fig. [S] The 
solid line is the negative, scaled probe force, amplified by a 
factor of 10; this force is weaker than the one in Fig.O (b) Re- 
sulting wave-packet dynamics in the Floquet representation, 
shown as a contour plot of \g\{k,t)\'^ . Under the first action 
of the probe force the wave packet is shifted in accordance 
with the generalized acceleration theorem p2|) . but not as far 
as in Fig. O such that it undergoes only partial Zener tran- 
sitions. After the dressing force is switched off, the second, 
reversed action of the probe force shifts the remaining part of 
the packet back to the Brillouin zone center, (c) Comparison 
of the initial wave packet (dashed line) with the part of the 
wave function that remains in the lowest energy band at the 
end of the process (solid line). 

ation theorem derived in Sec. IIIII takes the same form as 
its historic antecessor formulated by Bloch^ but applies 
to single quasienergy band dynamics, which can be dras- 
tically different from single energy band behavior. There 



are further features which can be carried over from the 
crystal-momentum representation to the Floquet picture 
and acquire a modified meaning there, such as the ex- 
pression for the group velocity of a wave packet or Zener 
transitions among different bands. 

The super Bloch oscillations considered in Sec. lIVI pro- 
vide a mainly pedagogical example which can be worked 
out in full detail analytically. Here the Floquet picture 
cannot exert its full strength, because one assumes a pri- 
ori that the driving force does not induce transitions 
from the initially occupied energy band to other ones, so 
that the historic acceleration theorem remains capable 
of describing the entire dynamics. The Floquet approach 
leads to exactly the same result, but implies a different 
viewpoint, separating the dc component of the force into 
one part which is resonant with the ac component, and 
together with the latter dresses the lattice, creating a 
quasienergy band; the remaining residual part of the dc 
force then probes this new quasienergy band, rather than 
the original unperturbed energy band. 

This theme of "dressing and probing" also prompts 
far-reaching strategies for achieving coherently controlled 
interband population transfer and even more. Two ba- 
sic examples for this have been discussed in Sec. |Vl but 
the possibilities obviously extend much farther. Utiliz- 
ing the generalized acceleration theorem, an initial wave 
packet may by split coherently into two components at 
an avoided quasienergy band crossing in a dressed lattice, 
and the lattice may then be redressed (that is, exposed 
to an ac force with different parameters) such that an- 
other quasienergy band structure is generated, possibly 
involving avoided crossings which affect only one of the 
daughter wave packets created in the first step, but not 
the other. Moreover, daughter wave packets can be made 
to move, possibly into different directions, and to inter- 
fere with other wavelets having been manipulated sepa- 
rately before in distant parts of the lattice. This vision 
apparently will be hard to realize with traditional solids, 
but it has come into immediate reach in current labora- 
tory experiments with weakly interacting Bose-Einstein 
condensates in driven optical lattices. Seen against this 
background, the generalized acceleration theorem almost 
provides a blueprint for a wave-packet processor. 



Acknowledgments 

M.H. thanks T. Monteiro for a thorough discussion 
of super Bloch oscillations. This work was supported 
by the Deutsche Forschungsgemeinschaft under Grant 
No. HO 1771/6. 



1 F. Bloch, Z. Phys. 52, 555 (1929). (Harcourt, Fort Worth, 1976), Chap. 12. 

^ N. W. Ashcroft and N. D. Mermin, Solid State Physics ^ C. Kittel, Quantum Theory of Solids, 2nd ed. (Wiley, New 



12 



York, 1987). 

* H. G. Roskos, M. C. Nuss, J. Shah, K. Leo, D. A. B. 
Miller, A. M. Fox, S. Schmitt-Rink, and K. Kohler, Phys. 
Rev. Lett. 68, 2216 (1992). 

^ J. Feldmann, K. Leo, J. Shah, D. A. B. Miller, J. E. Cun- 
ningham, T. Meier, G. von Plessen, A. Schulze, P. Thomas, 
and S. Schmitt-Rink, Phys. Rev. B 46, 7252 (1992). 

® J. Feldmann, Block Oscillations in a Semiconductor Su- 
perlattice, in Festkorperprobleme, Advances in Solid State 
Physics Vol. 32, edited by U. Rossler (Vieweg, Braun- 
schweig/Wiesbaden, 1992), p. 81. 

C. Waschke, H. G. Roskos, R. Schwedler, K. Leo, H. Kurz, 
and K. Kohler, Phys. Rev. Lett. 70, 3319 (1993). 

* A. A. Ignatov, E. Schomburg, K. F. Renk, W. Schatz, J. 
F. Palmier, and F. MoUot, Ann. Physik 3, 137 (1994). 

^ J. Rotvig, A.-P. Jauho, and H. Smith, Phys. Rev. Lett. 74, 
1831 (1995). 

A. A. Ignatov, E. Schomburg, J. Grenzer, K. F. Renk, and 
E. P. Dodin, Z. Phys. B 98, 187 (1995). 
" A. M. Bouchard and M. Luban, Phys. Rev. B 52, 5105 
(1995). 

^2 K. Unterrainer, B. J. Keay, M. C. Wanke, S. J. Allen, D. 
Leonard, G. Medeiros-Ribeiro, U. Bhattacharya, and M. 
J. W. Rodwell, Phys. Rev. Lett. 76, 2973 (1996). 
V. G. Lyssenko, G. Valusis, F. Loser, T. Hasche, K. Leo, 
M. M. Dignam, and K. Kohler, Phys. Rev. Lett. 79, 301 
(1997). 

^'^ G. Lenz, L Talanina, and C. M. de Sterke, Phys. Rev. Lett. 
83, 963 (1999). 

M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. 
Salomon, Phys. Rev. Lett. 76, 4508 (1996). 
^® O. Morsch, J. H. MuUer, M. Cristiani, D. Ciampini, and 
E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001). 
H. Jones and C. Zener, Proc. R. Soc. Lond. A 144, 101 
(1934). 

^® C. Zener, Proc. R. Soc. Lond. A 145, 523 (1934). 

J. Callaway, Quantum Theory of the Solid State (Academic 

Press, New York, 1976). 
^° M. Holthaus, Z. Phys. B 89, 251 (1992). 

21 K. Drese and M. Holthaus, Chem. Phys. 217, 201 (1997). 

22 S. Arlinghaus and M. Holthaus, Phys. Rev. A 81, 063612 



(2010). 

J. Zak, Phys. Rev. Lett. 71, 2623 (1993). 

V. Grecchi and A. Sacchetti, Phys. Rev. B 63, 212303 

(2001). 

25 B. Simon, Phys. Rev. Lett. 51, 2167 (1983). 
2'^ M. V. Berry, Proc. R. Soc. Lond. A 392, 45 (1984). 
2^ H. Sambe, Phys. Rev. A 7, 2203 (1973). 
2® E. Haller, R. Hart, M. J. Mark, J. G. Danzl, L. Reich- 
soUner, and H.-C. Nagerl, Phys. Rev. Lett. 104, 200403 

(2010) . 

28 A. R. Kolovsky and H. J. Korsch, J. Sib. Fed. Univ.: Math. 

Phys. 3, 311 (2010); arXiv:0912.2587 
^° A. Alberti, V. V. Ivanov, G. M. Tino, and G. Ferrari, Nat. 

Phys. 5, 547 (2009). 
^1 K. Kudo and T. S. Monteiro, Phys. Rev. A 83, 053627 

(2011) . 

^2 W. Kohn, Phys. Rev. 115, 809 (1959). 
W. V. Houston, Phys. Rev. 57, 184 (1940). 
A. Eckardt, M. Holthaus, H. Lignier, A. Zenesini, D. 
Ciampini, O. Morsch, and E. Arimondo, Phys. Rev. A 79, 
013611 (2009). 

M. Grifoni and P. Hanggi, Phys. Rep. 304, 229 (1998). 
^® S. Longhi, Phys. Rev. A 71, 065801 (2005). 

K. Kudo, T. Boness, and T. S. Monteiro, Phys. Rev. A 80, 
063409 (2009). 

^® A. Zenesini, H. Lignier, D. Ciampini, O. Morsch, and E. 

Arimondo, Phys. Rev. Lett. 102, 100403 (2009). 
^8 A. Eckardt, P. Hauke, P. Soltan-Panahi, C. Becker, K. 

Sengstock, and M. Lewenstein, EPL 89, 10010 (2010). 
'^^ J. Struck, C. Olschlager, R. Le Targat, P. Soltan-Panahi, 

A. Eckardt, M. Lewenstein, P. Windpassinger, and K. 

Sengstock, Science Express, published online 21 July 2011 

[D OLIO. 1126 /science. 1207239] 
■^^ O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 

(2006) . 

'*2 L Bloch, J. Dahbard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

■'^ A. Fratalocchi and G. Assanto, Phys. Rev. A 75, 013626 

(2007) . 

''^ S. Arlinghaus and M. Holthaus (unpublished). 



